Neutrino emission due to Cooper-pair recombination in neutron stars revisited 
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I. INTRODUCTION 

In minutes/hours after the birth a neutron star cools 
down to a temperature T ~ MeV via a neutrino transport 
to the surface and then becomes transparent for neutri- 
nos. Hereafter during ~ 10 5 yr the cooling is determined 
by emissivity of neutrinos produced in direct reactions 
[E H, S 0) IE El- For such temperatures neutrons and 
protons in the neutron star interior are highly degener- 
ate. Therefore the rate of the neutrino production is 
suppressed by the reaction phase space, the stronger, the 
more nucleons are involved. 

The most efficient are one-nucleon processes, e.g., n — > 
pep, called direct Urea (DU) reactions. Their emissivity 
is e DU ~ 10 27 x T$ (n/n p 3 8(n - n° u ) JSfc, see [f, 
where Tg = T/(10 9 K) , n is the nucleon density measured 
in units of the nuclear matter saturation density no- The 
DU processes are operative only when the proton fraction 
exceeds a critical value of 11-14%. Equations of state 
constructed from realistic nucleon-nucleon interactions, 
like Urbana-Argonne one [8j, show that this condition is 
fulfilled only at very high densities. This implies that 
the DU processes may occur only in most heavy neutron 
stars, e.g., with masses ~ 2M for the equation of state 
Q, where M Q = 2 x 10 33 g is the solar mass. At n ~ n 
the proton fraction is typically about 3-5%, cf. Fig. 2 in 
Ref. [|. 

In the absence of the DU processes, most efficient be- 
come two- nucleon reactions, e.g., nn — > npeD, called 
modified Urea processes (MU) with the emissivity £ MU ~ 
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10 21 x T 9 8 (n/n ) 2/3 , cf. Note the smaller nu- 
merical prefactor and the higher power of the tempera- 
ture for the emissivity of the two-nucleon processes com- 
pared to the DU emissivity. In-medium change of the 
nucleon-nucleon interaction in the spin-isospin particle- 
hole channel due to pion softening may strongly increase 
the two-nucleon reaction rates, which, nevertheless, in 
all relevant cases remain significantly smaller than that 
for the DU 0,[ll|. The nucleon bremstrahlung reactions, 
like nn — > nnvv (nB) and pp — > ppvv (pB), have an order 
of magnitude smaller emissivity than MU. 

At low temperatures the nucleon matter is expected to 
undergo a phase transition into a state with paired nucle- 
ons [12j . The neutron superfluidity and/or proton super- 
conductivity take place below some critical temperatures 
T c .„ and T c>v , respectively, which depend on the density. 
At densities n < (2-4) no neutrons are paired in the ISo 
state and in the 3P2 state at higher densities. Protons 
are paired in ISo state for densities n < (2-4)no- Paring 
gaps, A i? are typically ~ (0.1-1) MeV and depend cru- 
cially on details of the interaction in the particle-particle 
channel, see Fig. 5 in 

The gap in the energy spectrum significantly reduces 
the phase space of the nucleon processes roughly by the 
factor exp(— A/T) for the one-nucleon DU process and 
exp(— 2A/T) for two-nucleon processes. However, even 
with inclusion of the nucleon pairing effects the DU rate 
is large enough that the occurrence of these processes 
would lead to an unacceptably fast cooling of a neutron 
star in disa gre ement with modern observational soft X- 
ray data (9l.Tl3l. [l4| . This statement has been tested with 
gaps varying in a broad band allowed by different mi- 
croscopic calculations. Certainly, DU processes could be 
less efficient if one kept gaps finite also at high densities. 
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Microscopic calculations do not support this possibility. 
Thus, according to the recent analysis the DU processes 
most probably will not occur in typical neutron stars with 
masses in the range of 1.0-1.5 Mq, based on the cooling 
and the population syntheses scenarios [H, HH EH] ■ 

The superfluidity allows for a new mechanism of neu- 
trino production associated with Cooper pair breaking 
or formation (PBF), e.g., the reaction n — ► nv v and 
p — > pv where one of the nucleons is paired. For ISo 
neutron pairing the PBF emissivity was evaluated first 



in 16 1 in the Bogolyubov ^-operator technique and then 
IZjEH within the Fermi-liquid approach. 



The proton PBF emissivity was estimated in [17j with 
account for in-medium renormalization of the nucleon 
weak-interaction vertex due to strong interactions. Mix- 
ing of electromagnetic and weak interactions through the 
electron-electron-hole loop can additionally change the 
proton vertex nam. There are also relativistic cor- 
rections to the axial-vector coupling vertices of the or- 
der v F t , vf i are neutron and proton Fermi velocities, 
i = n,p, .21] . These three effects together resulted in 
a one or two orders of magnitude enhancement of the 
proton PBF emissivity compared to that evaluated with 
the free vector-current vertex. Thus, one concludes that 
neutron PBF and proton PBF emissivities can be equally 
important for neutron star cooling depending on the pa- 
rameter choice and, especially, on relation between gaps 
A p and A„. The PBF emissivity for the 3P2 neutron 
pairing has been analyzed in [22I ]. 

Following 0, [III IH1 the emissivity of the neutron PBF 

and proton PBF processes is estimated as e lPFB ~ 10 28 x 
(A l /MeV) 7 (T/A. i ) 1 / 2 (n l /n ) 1 /3 e -2A I /T for T « 

Aj and i = n,p. Having a large numerical pre- factor 
and very moderate temperature dependence of the pre- 
cxponent, these reactions significantly contribute to the 
neutron star cooling provided gaps are not too small. 
These processes have been included in the cooling code 
rather recently [23]. From that time the PBF reactions 
are the main part of any cooling-scenario together with 
the MU processes [H, Q, [H, [25[ . Uncertainties in the 
pairing gaps are large. Therefore surface temperatures 
of neutron stars computed in different approaches vary 
significantly. 

Kundu and Reddy [2(| and Leinson and Perez [27| 
made an important observation: all previous calculations 
of the neutrino reactions in superfluid matter disrespect 
the Ward identity and, as a consequence, the conserva- 
tion of the electro- weak vector current. 

The Ward or in general case Ward-Takahashi identities 
impose non-trivial relations between vertex functions and 
Green functions, which synchronize any modification of 
the Green-function with a corresponding change in the 
vertex function. Satisfying these relations one assures 
that the symmetry properties of the initial theory are 
preserved in actual calculations. For instance we start 
with the theory of weak interactions with a conserved 
vector current. The current would remain trivially con- 
served in calculations with only bare vertices and bare 



Green functions. In strongly interacting systems the 
Green functions change necessarily, but for quasiparti- 
cle Green functions the current conservation is easily re- 
stored by a proper inclusion of short-range correlations 
in the vertices, cf. development of the Fermi-liquid the- 
ory by Migdal flEl [28l [29| . Following these two simplest 
cases, Refs. [l6l I2l|j22| did not incorporate any medium 
effects, whereas |17|used dressed quasiparticle Green 
functions together with dressed normal vertices. In case 
of the superfluid system situation is more peculiar. Since 
in the superfluid system the nucleon Green function no- 
toriously differs from the free one, the vector-current ver- 
tex must get corrections even if no other interaction be- 
tween quasiparticles is included. Additional anomalous 
vector-current vertices disregarded in previous calcula- 
tions must be properly accounted for. These corrections 
cancel exactly the vector-current contributions to the neu- 
trino emissivity for zero neutrino momenta, cf. (27| . 

Assuming that the axial-vector current contributes 
only little to the PBF emissivity, Ref. 
the PBF emissivity calculated in [lj| 



27 
17 



J20 ~ 10- 



claimed that 
_, Sl| is to be 
3 for n ~ rig for 



suppressed by factor ~ 
neutrons and by ~ 10~ 7 for protons in the case of ISo 
pairing. Such a severe reduction of the neutron and pro- 
ton PBF emissivities could significantly affect previous 
results on the neutron star cooling dynamics. Ref. [3lj] 
revises results [27j applying expansion in q 2 parameter 
and putting vp.n = 0. Ref. [3 II ] claims that suppression 
factor for the neutron PBF emissivity is ~ T/m*, where 
m* is the nucleon effective mass. This would reduce the 
neutron PBF emissivity by the factor ~ 5 x 10~ 3 for 
temperatures T ~ 0.5T Cj „, cf. Fig. 5 in [3l| . 

Refs. [13, [3l| used the convenient Nambu-Gorkov ma- 
trix formalism developed to describe metallic supercon- 
ductors [H, HH. The price paid for the convenience is 
that the formalism does not distinguish interactions in 
the particle-particle and particle-hole channels. Such 
an approach is, generally speaking, not applicable to 
the strongly interacting matter present in neutron stars. 
In the nucleon matter at low temperatures the nn and 
pp nucleon-nucleon interactions in the particle-particle 
channel are attractive, whereas in relevant particle-hole 
channels they are repulsive [ij], [H, H(|. The adequate 
formalism was developed by Larkin and Migdal for Fermi 
liquids with pairing at T = in Re f. [ 341) and generalized 
then by Leggett for T ^ in Ref. (35^ 

In the present paper using Larkin-Migdal-Leggett for- 
malism we analytically calculate neutrino emissivity from 
the superfluid neutron star matter with the ISo neutron- 
neutron and proton-proton pairing. Both normal and 
anomalous vertex corrections are included. We explicitly 
demonstrate that the Fermi-liquid renormalization [341 ] 
respects the Ward identity and vector current conserva- 
tion. Our final estimations of neutron and proton PBF 
emissivities differ from those in [53, Hl| . We find that the 
main term in the emissivity ~ v F i follows from the axial- 
vector current, whereas the leading term in the emissivity 
from the vector current appears only at the v F i order, as 
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in |22| . 

In next section we present the general expression for 
the emissivity of the neutron PBF processes formulated 
in terms of the imaginary part of the current-current 
correlator for weak processes on neutral currents. Then 
within the Fermi-liquid approach to a superfluid we in- 
troduce Green functions and the gap equation. In Sect. 3 
we formulate and solve Larkin-Migdal equations for ver- 
tices and apply them to calculate the imaginary part of 
the current-current correlator. In Sect. 4 we calculate 
the neutrino emissivity of the neutron PBF processes. 
Also in this section and in Appendix a comparison of our 
results with results of previous works is performed. In 
Sect. 5 we consider emissivity of the proton PBF pro- 
cesses. Conclusions are formulated in Sect. 6. 



II. GENERAL EXPRESSIONS: EMISSIVITY, 
GREEN FUNCTIONS AND PAIRING GAPS 



A. Neutrino emissivity 

The weak neutrino-neutron and neutrino-proton inter- 
actions on neutral currents are described by the effective 
low-energy Lagrangian 



C = 



G 
2\/2 



(1) 



where 1^ = v^j^ (1 — 75) v is the lepton current and = 

gy *i 7 P and A? = 7^75 stand for nucleon 

(neutron or proton) vector and axial- vector currents with 

nucleon bi-spinors The coupling constants are gy = 
gv = -1, g W = c v = 1 - 4sin 2 6> w ^ 0.04 and g^ = 
—g^ = Qa = 1-26. The Fermi constant is G « 1.2 x 
10~ 5 GeV~ 2 . For the non-relativistic nucleons w 

^l(p') (1. (P + P)/2™) MP) and A» « Vl(p') (°(p + 
p)/2m, a) ipi(p) j where a = (ci, 02, 03) are the Pauli 
matrices acting on nucleon spinors ipi > ancl P an d P are 
outgoing and incoming momenta, m is the mass of the 
free nucleon, cf. [36| . 

Neutrino emissivity for one neutrino species can be cal- 
culated as 



G 2 



d3qi ^ ,/b(,)2^ X ( 9 ) 1 (2) 



(27t) 3 2lji (27r) 3 2w 2 



where q = (w,<f) = qi + q 2 , qi.2 = (wi, 2, 91,2) 
are 4-momenta of outgoing neutrino and antineu- 
trino, /b(w) = l/(exp(u;/T) — 1) are Bose occupa- 
tions, and 3% is the imaginary part of the suscepti- 
bility of the nucleon matter to weak interactions, i.e. 
the Fourier-transform of the current-current correlator 
(faWHx) - A^x)l»{x))(V»{y)ll{y) - A" (y)ll(y))) , 
for weak processes. The sum in @ is taken over the 
lepton spins. 



According to the optical theorem 3% can be expressed 
as the sum squared matrix elements of all available re- 
actions with all possible intermediate states, \M\ 2 . A 
particular contribution to ^ \M\ 2 can be also calculated 
within the Bogolubov ^-operator approach for a given 
form of the nucleon-nucleon interaction, as it has been 
done in Refs. [22J. In this approach, however, an 
account of further in-medium modifications of nucleon 
propagators and interaction vertices is obscured by a 
danger of double counting. The Green-function tech- 
nique for Fermi liquid (2&I30I] is more suitable for such 
extensions as it was demonstrated in [Til [l7j . We will 
follow the Green-function approach for superfluid Fermi 
liquids 0, [34|]. As a simplification we will focus on 
the low temperature limit T <C A. The temperature 
dependence enters through nucleon occupation factors 
tx e" A / T (l + 0(T 2 /A 2 ) and also as 1 + 0(T 2 /e|) cor- 
rections in the low-temperature expansion of standard 
Fermi integrals, when the high energy region e > A is 
dominating in the integrals. Since the boson occupation 
factor fs in (J2J> generates already the leading exponent 
e -2A/T we can eva i ua t e for T = 0, see also discussion 
below. 



B. Nucleon Green functions and pairing gaps 

The nucleon Green function for the interacting system 
in a normal state ("n.s."), i.e. without paring, is given by 
the Schwinger-Dyson equation, which in the momentum 
representation reads 

Gn.sip) = G (p) + Gob) %,...(?) G n . s .(p) 

with Go (p) — Go (e, p) 1 = 1/ (e — e p + i sng e) , where 
1 is the unity matrix in the spin space. All informa- 
tion on the interaction is incorporated in the nucleon 
self-energy S n . s ., being a functional of the Green func- 
tions G n . s . . In absence of the spin-orbit interaction the 
full Green function is also diagonal in the spin space, i.e. 
G n . s . = G n . s . 1 • For strongly interacting systems, like 
dense nucleon matter, the exact calculation of G n . s . is 
extremely difficult task. However, for strongly degener- 
ate nucleon systems at temperatures, T, much less than 
the neutron and proton Fermi energies, eF.i, i — Ti,p, 
fermions are only slightly excited above the Fermi see. 
So, the full Green function of the normal-state is given 
by the sum of the pole term and a regular part: 



G„. s .(p) = 



i sgne 



+ G reg(p) , 



(3) 



where the excitation energy is counted from the nu- 
cleon chemical potential /1, e p = p 2 /(2m*) — jj, , /u, ~ 
e F = Pp/C^m*) for low temperatures under considera- 
tion, p-p is the Fermi momentum. The effective mass and 
the non-trivial pole residue are determined by the real 
part of the self-energy, as a -1 = 1 — (9KS n . s ./9e) F and 
1/m* = a (l/m + 2d i RJ: n . s ./dp 2 ) w . The subscript "F" 
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indicates that the corresponding quantities are evaluated 
at the Fermi surface (e, e p — > 0). According to Ref. [29[ 
only the pole part of G„. s . is relevant for the description of 
processes happening in a weakly excited Fermi system. 
The regular part can be absorbed by the renormaliza- 
tion of the particle-particle and particle-hole interactions 
at the Fermi surface. The quantities m* and a can be 
expressed through the Landau-Migdal parameters char- 
acterizing the fermion interaction at the Fermi surface at 
zero energy-momentum transfer. The imaginary part of 
the self-energy, QE n . s ., can be omitted in the pole term 
of the Green function ^ in the low-temperature limit 
(quasiparticle approximation) . 

In a system with pairing a new kind of processes such 
as transition of a particle into a hole and a condensate 
pair and vice versa become possible. The one-particle 
one-hole irreducible amplitudes of such processes can be 
depicted as [34[ 



;A( 2 ) 



(4) 



Besides the "normal" Green functions shown by the 
" thick" line iG = — » — one introduces " anomalous" 
Green functions iF^ = -» — «- , and iF^ = — « — *— . 

The full normal and anomalous Green functions are re- 
lated by Gor'kov equations 

G(p) = G n . s .( P ) + G n .J P ) A«(p) F™(p) , 
F^(p) = G h n . s .(p)A^(p)G(p). (5) 

The second equation involves the normal-state Green 
function of the hole (superscript "h"), which in the ab- 
sence of a spin-orbit interaction is simply i Gjj s (p) = 
i G n . s .(—p) = —* — • In the case of the ISq pairing the 
spin structures of the anomalous Green functions and 
the transition amplitudes ^ are simple: A^ 1 ) = A^ 2 ) = 
A i cr 2 and F^ = F^ = Fia 2 . Eqs. © are to be com- 
pleted by the equation for the amplitude AW(p), 

^ =/^^^^')]: d C [G(p')A(V)GL.(p')]e>(6) 



energies far off the Fermi surface. One may renormal- 
ize [29|, [3(| the interaction (V — > T^) in such a manner 
that integrations go over the region near the Fermi sur- 
face and only the quasiparticle (pole) term in the Green 
function is operative. Advantage of the Fermi-liquid 
approach is that all expressions enter renormalised ampli- 
tudes rather than the bare potentials. For \p\ ~ pp — \p | 
the effective interaction amplitude is a function of the an- 
gle between p and p only. The amplitude in the particle- 
particle channel is parameterized as 

[f «] 2 = 4 (n, ff) (i<7 2 )?(i<7 2 )5 + rf {ft, n') (ia 2 a) a b {a ia 2 ) d: 
and the interaction in the particle-hole channel is 

[fX c = r o (*> % d d + r i ft ^ ft)S- 

Here and below ft = p/\p\ and n' — p' /\p' \ . Superscript 
V indicates that the amplitude is taken for \qvp \ <C u 
and u> -C ep, where lu and q are transferred energy and 
momentum. Amplitudes Tq, are expanded in the 
Legendre polynomials. 

Integrating over the internal momenta in loops we 
can separate the part accumulated in the vicinity of the 
Fermi surface / ~ f x f d<S> p with Jd$ p = 

P = ^ being the density of 
states at the Fermi surface. After the Fermi-liquid renor- 
malization ^ reduces to 

A(n) = -AoQ* (H, ft) (8) 
A = f d$ p G n . s .(p)Gi s .(p)8(t-e p )na 2 P H2Z/A), 

where we denoted (. . = J -r-r (■ ■ ■) and £ ~ £f • One 
usually determines the gap supposing £ = ep . 



III. CURRENT-CURRENT CORRELATOR, 
EQUATIONS FOR VERTICES, AND VECTOR 
CURRENT CONSERVATION 

A. Current-current correlator 



where V stands for a two-particle irreducible potential, 
which determines the full in-medium particle-particle 
scattering amplitude. The potential V can be separated 
in the scalar and spin-spin interactions defined as 

[V]Z = V v {ia 2 y b (ia 2 ) d + V A {ia 2 a)t {ffio 2 ) d . 

Solution of the Gor'kov equations is straightforward. 
The relevant pole parts of Green functions are 



G(p) = £ 2 -£ 2 + z0sgn £ ' f(p) 



-a A 



e 2 — _E 2 + iOsgne 



, (7) 



where i? 2 = e 2 + A 2 . Integrations over the internal 
momenta in fermion loops, e.g., over p' in ©, involve 



Applying theory of Fermi liquids with pairing [29l . |34| 
we can present contributions to the susceptibility \ i n 
terms of the diagrams 



IX = 



0-0*0*0- 



.(9) 



Here dash line relates to the Z-boson coupled to the neu- 
tral lepton currents, vertices on the left are the bare ver- 
tices following from the Lagrangian (JTJ) . The right-hand- 
side vertices f, f h , f^ 1 ' and are the full vertices de- 
termined by the diagrams shown in Fig. [T] The blocks in 
Fig. [T] correspond to the two-particle irreducible interac- 
tion in the particle-particle channel, T^, and the particle- 
hole irreducible interaction in the particle-hole channel, 




FIG. 1: Graphical representation of eqs. (|13[) 



T w . We emphasize that only chains of bubble diagrams 
are summed up in this particular formulation. Thus, the 
imaginary part of x accounts only for one-nucleon pro- 
cesses. In order to include two-nucleon processes within 
a quasi-particle approximation one should add diagrams 
with self-energy insertions to the Green functions and 
iterate the Landau-Mi gda l amplitudes T^'t in Fig. [1] in 
the horizontal channel [l7l . |37| . In general case for parti- 
cles with widths the interpretation of different processes 
contributing to 3% is more peculiar and needs another 
re-summation scheme (37j . 

Taking the imaginary part of x we cu t & diagram 
through two fermion lines. Cuts of neutron lines cor- 
respond to neutron PBF processes and those of proton 
lines correspond to proton PBF processes. Since the neu- 
tron density in a neutron star is much higher than the 
proton density, we can drop all diagrams where proton 
lines are uncut from the set of the bubble chains included 
in ©. 

Refs. [ill Hi], [22| considered only first two diagrams 
in © with bare vertices. Refs. @,Ej| treated those two 
diagrams with full vertices on the right, whereas one must 
consider all four diagrams with the full vertices. 

Vector and axial-vector currents contribute to x sepa- 
rately, i.e., x = Xv + XA, where 

+ G + r^F™ + F« tM6_},a = V,A. (10) 

Here and below we use the short-hand notations G± = 
G(p±q/2) and the analogous one for the F± Green func- 
tions. All left vertices in §§§ are the "bare" vertices, 
t% , after the Fermi-liquid renormalization [2^, [3(| , = 
[1 + Tq (G+G-Y] t°, which involve the particle-hole ef- 
fective interaction, Tq, integrated with off-pole parts of 

the Green functions (G+G_)" = lim^o / (Hylq G+G_, 

and r° follows from (fTJ). The difference between r° and 



t£ can be cast [29] in terms of a local charge of the quasi- 
particle e a — ar^ /r° . Then 

T v = 9v (t v o l Q - t v1 I ) , t v o = — , r v l = — v , (11) 

T A ~9A yA.l - T A fi <7l ), T A Q = — , T A1 = — V. 

For the vector current ey = 1 and the vertices Ty and 
satisfy the Ward identity wTy — qfy 1 = Gi P °. le ^ _1 (p + 
g/2) — Gn.s.' 1 (p—q/2), with the pole part of the normal 
state Green function, Gn.s!^ = G n . s . — G rog . The local 
charge for the axial-vector current differs from the unity 
varying in different parameterizations as &a — 0.8-0.95, 
as it follows from studies of the Gamow- Teller transitions 
in nuclei, see [29|, [H, [39| and references therein. 



B. Larkin-Migdal equations for full vertices 

Consider first one sort of nucleons, e.g., neutron. At 
the Fermi surface the full vertices f, f h , f W and f ( 2 ) can 
be treated as functions of out-going momentum q and the 
nucleon Fermi velocity v = vfh, n — p/p. Their general 
structures are 

T v = 9v (tv,o lo — tv,i l) j t v = gv (tv,qIq + Tv,il) , 

(1) (2) /- , ^ j\ . 

t v = -t v ' = -g v [Tv,olo ~ tva I ) l&2 , 

TA — — gA (ta,1^1q — TA,0 Vl) , 
T A = — gA\ — TA,W h—TAfiCF I), 
T A = +gA(TA,l&lo-TAfivi)iV2, 

t { a = -gAi(T 2 (ta,i&Io -TAfiSl) . (12) 

Superscript "T" denotes matrix transposition. 

As follows from the diagrammatic representation of 
Fig.[T]the full vertices obey Larkin-Migdal equations [34j : 



T a ,a{n,q) = r£ (n,g) + (T"(n,n') [L(ri ' , q; P a ,o) T afi (ri \q) + M(ri \q)T afi {ri \q)]) H , , 

r a ,o(n, q) = - <r«(n, n') [(N{n>, q) + A ) T a , {n', q) + 0(n>, q; P a<0 ) T a , (n', q)] ) n , , 
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T a ,i(n,q) = T%,i{n,q) + (r£(n,n') [L(n>, g; P a>1 ) f 0)1 (n', q) + M {n>, g) r a>1 {n', g)] > fl/ , 

T a ,i(Ti,g) = - <ri(n,n') [(JV(n', g) + A ) f«,i(n', g) + 0(n', g; P ,i) w(< «)] )# 

I 



(13) 



where a = V,A and Py i0 = — Py,i = — Pa.o = Pa,i 
I . In order to write the one set of equations for both 
vector and axial-vector weak currents we introduced a 
new notation for the effective interaction = Tq, if 
a = V, and = T^, if a = A. Functions L, M, N, 
and O arc defined as 



L(n, g; P) = /d$ p G+G_ - (G+G_) w - P+F_P 
qv 



a 2 p 



M(n,q) = /d$ 



N(n, q) 



uj — q v 
G+F_ - F + G- 



2 u + qv 

- a P rt A g( Z )> 



.(l- fl (z))- fl («)(l+P)/2 



2A 



d$, 



G + G_ - (G p G p )6(£ 



F>F_ 



2 - (g?) 2 
a P — 9{z), 



4 A 2 

0(n,q;P) = - /d$ p G+F_ + F+G h _P 
'uj + qv uj — qv 



a 2 p 



4A 



+1/2 



g{z 2 ) 



4A 



d;r 



P 



4 z 2 x 2 - z 2 + 1 + i 



(14) 



(15) 



-1/2 



arcsinhV z 2 — 1 in9(z 2 — 1) 



2z\/i 2 ^T 



2 ^ 2 - (gV) 2 
Z = 4A 2 ' U = VF 



Expressions (fl4|) . (|15|) are derived for T = 0. For fi- 
nite temperatures, T < T c , all expressions in (|14[) . ex- 
cept for L, hold as well, but with g(z 2 ) — > g(z 2 ,T) and 
A — > A(T). Generalization of expression for L requires 
introduction of one more temperature dependent inte- 
gral besides g. Such expressions were derived by Leggett 
in [35|. As follows from these expressions there arises 
an essential simplification in the limit of low tempera- 
tures, T <c A. We exploit the fact that to calculate 
the PBF emissivity we need only imaginary part of the 
current-current correlator 3%. Since uj > 2A for the PBF 
kinematics, the emissivity is exponentially suppressed by 
e -2A/T stemming from the Bose occupation factor /b(w) 
in ©. Therefore, we may take 3% K ^9{T = 0) since 
it is already multiplied by the term vanishing for T — > 0. 
Not accounted temperature corrections in 3% prove to be 
~ l + 0(e- A / T (l + T 2 /A 2 )) + 0(T 2 /e|). The latter term 
follows from the expansion of Fermi integrals when the 



integration goes over energy regions far from the Fermi 
surface. Such corrections are small in the limit T <C A 
and we omit them. 

Using vertices (|llll2p in (fT0|) the correlators \v and 
Xa can be expressed as 

Xv(q)=gy ({ l o -vl) {laXv,o{n,q) ~ Xv,i{n, q) I t )) fJ , 

XA(q)=9A((lov-l) (llxA,i{n,q) - XA,o(fi, q)l )) n , 
Xa,o(n, q)=L(ft, q; P a ,o) r 0) o(n, q) + M(n, q) r a>0 (n, q), 
Xa,i(n, q)=L(n, q; P 0)1 ) f ,i (ft, q) + Af (n, q) T ,i(n, q). 



Solution for vector and axial-vector parts of the 
current-current correlator 



It is natural to expect that first and higher Legen- 
dre harmonics of Y^'^in, n') are smaller than the zero- 
th ones due to the centrifugal factor [291 ] . This al- 
lows us to retain only zero harmonics T^'^(ft,Ti f ) = 



const , expressed through dimension-less Landau- 

f^/(a 2 p(n )) and 



Migdal parameters as [i 
r^'* = g u I (a 2 p(no)) ■ The values of parameters are 
extracted from the analysis of atomic nucleus experi- 
ments [29T |38|, |39[ or in some approximations can be 
calculated starting from a microscopic nucleon-nucleon 
interaction (40| . Actually, in isospin asymmetric mat- 
ter f u and g u are different for interactions between two 
neutrons (/£ n , g% n ), two protons (/£,, g% p ) and neutron 
and proton {f% p , gf lp ). Note that values f% n , f£ p are 
necessarily positive, the requirement of the stability of 
the nucleon matter, whereas corresponding values in the 
particle-particle channel /|„, /| p are negative, otherwise 
there would be no ISq pairing. In this respect our deriva- 
tions differ from those which do not distinguish interac- 
tions in particle-hole and particle-particle channels and 
use Nambu-Gorkov formulations with one bare potential 
(V < in our case). 

For the angular-independent amplitudes (only zero- 
th harmonics are included) the Larkin-Migdal equations 
(ITcH) get simple solutions: 

7-0,0 (g) = ia(q; Pafi) T% fi , 
7 - 1 (g;P) = l-r- (C(n,q;P)) H , 



C(n,q;P)=L(n,q;P) 



(Q(n 7 q;P)) n 
(N(n,q)) n 



M(n, q), 



1~a,0 



(<?) 



(Q(n,q;Pg.o))n 
(N(n, q)) n 



r a,o(q) ■ 



(16) 



We have exploited here the relation 1 = — Tq (A ) fol- 
lowing from the gap equation {5}. Although integrals 
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Xa,i(n,q) = vj a (q;P atl )jC(fi,q;Pa,i) + 5Xa,i(n, q) , 
<%,i(n, q) 



M(n,q) ,-. 

■ {0(ri ,q\P a ,i){v - v ))„-, 



+C(n, q; P a ,Mq; P a<l ) (C(n',q; P a , x ){v' - v))-, . 

Here j a are precisely those nucleon-nucleon correlation 
factors that have been introduced in (l7| . They depend 
on Landau-Migdal parameters in the particle-hole reac- 
tion channels. 



D. Vector current conservation 

Now we are in the position to verify that the correlator 
of the vector current Xy supports the current conserva- 
tion. First we note that there are convenient relations 

(C(n,q;P)-C(n,q;P)) R = 0, 
(uC(n,q;±l) - C(n,q;Tl) (qv)) a = 0, 
(qXa,i(n, q))n = -fa(Pa,i)((qv) £(n, q; P a .l))n, 
((qv) (u£(n,q;+l) - (qv ) C(n, q; -1) 
-q5xv,i(n,q,^ = 0))) H = q 2 a 2 n/m* . (19) 

These relations help us to establish important properties 
of the vector current correlators (I18P : 



(wxv,o -qxvs)n = 7v{q;+i)iv(q;-i)uTy 

x (C(ft t q; +1)) ((C(n, q- -1)) - (jC(n, q; +1))) 
= 0(rgq 6 v^/u J 6 ), 
%((qv) (uxv,o -qXv,i))a = wTy (C(H,q;+l)) 
x ((qv ) [w{q\ +l)£(n, q; +1) - j v (q; -l)C(n, q; -1)]) 

= 0(rgfv%/w 6 ). (20) 

We use here expansion of C and C in the series for 
\q\ Vf/uj <C 1: 

C(n,q;+1) = -J^— + g y 2 (l - x 2 ) 
1 — y x 

x (l + yx + y 2 {\+x 2 ))+0{gy b ) 1 

C(n,q;-1) = — — gyx(l + yx + y 2 x 2 + y 3 x 3 ) 

1 — y x 



C(n,q;-1) = 



in (fT"3|) do not produce terms oc v for constant T%'£, + 
the vector vertices f a .i and T a ,i gain new terms propor- ^. = 

tional to q, thus, f a ,i(n, q) — r^ ^n, g) + n q T^\(q) and 
T aA (n,q) = n q T^l(q), where n q = q/\q \ and 

T^l(q)=j a (q;Pa,l)K(C(fi,q;P a>1 )(nn q ))n, (17) 

^{l)( n \ - (0(H,q;P a ,l)) Jl) _ (Q(H,q;Pa,l)(nn q )) H 
T a,l\q) ~ {N(n,q)) H T a,l (N(n,q)) H ' 

C(n, q; P) = L(n, q; P) - 0(n, q; P) . 

With (I11I16I17|) we cast X >Z = (Xa,o,Xa,i) as 

Xa,o(n, q) = -fa(q\ Pafi) C(H, q] P a fi) , (18) 



0(gy 5 ). 
yx 

i — y x 

9V 2 (\~ 
yx 



gyx{l + y 2 x 2 ) 



x 2 )(l 

- gy 2 x 



y 



2 a+x 2 )) 



0(gy 5 



2 ~ 2 (l + y 2 x 2 ) 



1-yx 
gxy 3 (l-x 2 ) + 0(gy 6 )., 
qv F /uj,x= (qv)/\q\ v F , 

2 



(21) 



fJ 



>'«')) 



Relations (j2"0|) demonstrate that the imaginary part of 
the vector current correlator calculated with full vertices 
(|16ll7p is transverse, 

^xvhqu = 0(rgq 6 v e F /u; 6 ), (22) 

at least up to terms of the higher order than 
/ w (|q| vp/uj) 5 g , which are beyond the Fermi-liquid ap- 
proximation for the Green functions (|3l7p . This ensures 
conservation of the vector current. Note that for r w = 
or g = we have 3(t^ Xy)n<iV = and then the vector 
current is conserved exactly. Since g oc A(T) and A = 
for T > T c , we have proven in passing that the vector 
current is conserved exactly above T c . 

In order to prove the transversality of the real part 
of the vector-current correlator it would be necessary to 

include the tadpole diagram contribution (^) , where 

the coupling originates from the " gauging" of the kinetic 

-2 

term, ip V ip/2 m* , of an effective non-relativistic nucleon 
Hamiltonian. 

Some comments on the approximations done in pre- 
vious works would be here appropriate. In all previous 
works the vector current contribution has been consid- 
ered as the dominant term for the case of ISq pairing. 
Expressions for the PBF emissivity in [l(| [U [UJ can 
be recovered, if we put Svy p(n, q) = 3L(n, q; +1) and 
i = 0. Result of @, LLD] is obtained by taking 
9Xv,ofc?) = ZL(n,q;+l)/(l + T%$tL(n,q;+l)) 2 and 
also ^Xv,i — ■ Setting = in the limit u> <C A and 
q = we reproduce expressions of Ref. [13] • 

Note that relations (p~9j) do not hold with C replaced by 
L and, hence, the transversality relation (|2"2")) is spoiled, 
if one ignors the anomalous vertex terms. 



NEUTRINO EMISSION VIA NEUTRON 
PBF 



IV. 



After correlators (fT8| are established it remains to take 
the sum over the lepton spins and integrate over the lep- 
tonic phase space in @. The latter can be easily done 
with the help of the Lenard integral [4l[ 

A2_ ^ { p^ }(5 (4) ((?i 



(2tt) 3 2o;i (2tt) 3 2u; 2 



?2 



48 7T 5 



q 



■q) 

2 ).(23) 
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Now the neutrino emissivity ([2]) can be cast as 



8 J J 6?r 4 

_* 2 few 
" 9a 



Qa(w) 



240 7T 4 y 

' d\q\q' 



d 3 qi d 3 q 2 
2u>i 2lu2 



dww"/B(w)Qa(w), (24) 



<5 (4) (gi + g 2 -g) 



(25) 



(26) 



In K a the sum is taken over the lepton spins. Shortening 
notations we introduced in ([24]) effective couplings g* a = 
&a 9a- The particular normalization of the quantity Q a is 
chosen so that for Q a (u>) — Q°{uj) with 

QW(«) - -p%(c 2 /4A 2 ) (27) 
we obtain expression for the neutron PBF emissivity 



oc 

e - ) = ^0^^ 1{z) = J 



<iyy 5 



-2zy 



,(28) 



which coincides with the old result [la . Il7| after the 
replacement e~ 2zy — > ^ eZ J~ +1 yz ■ From now on we restore, 
where it is necessary, subscripts "n" or "p" to distinguish 
neutron and proton PBF processes, respectively. 



A. Emissivity on vector current 

For the vector current we have 

Q 2 (Xv,o(n,q))n + ((qv) qxv,i(n,q))n 



+ {u 2 ~ q 2 ) (vxv,i(n,q))n 
- ^ ((qv)xv,o(n,q))n- (qxv,i(n,q))n -(29) 



Using relations (118)) and ([20]) we can simplify ([29]) as 

k v = (q 2 - cj 2 ) 3?(xko(^,<?) ~ vxv,i(n, q))n- (30) 

Both scalar and vector components in (|30[) are of the 
order Vp, 



%(xv,o(n, q)) 
^(vXv,i(n, q)) 



a — *4 4 

4g«| 2 

45^ 4 a P ^UA2 



2 



2gV 
9cl , 2 - ^~^ 4 A 2 



We have put 7v — > 1 since — 1 + 0(f^ n vp). 



Note that the first term in ([30]) . cx 3xv,0: would give a 
negative contribution to Qv- Only due to the presence of 
the vector component of the vertex, second term in ([30]) , 
the full expression for the reaction probability becomes 
positive. This is because we used Ward identities, which 
impose relations between zero- and vector components. 
However, if one keeps in ([2"9"| only the first term related 
to the zero-th component of the vertex and drops other 
terms, as it was done in early works, the expression for 
the reaction probability would be also positive. 

Then in terms of Qv we get 



(31) 



Finally for the neutron PBF emissivity on the vector cur- 
rent we obtain (for one neutrino flavor) 



nPFB 

t vv,V 



(On) 2 J* 4 
t vv 9V gi V F,n , 



(32) 



Note that in spite of Ref. [27[ used approximation uj <C 
A„, which is not fulfilled in PBF case, our expression 
([3"2"]) only slightly deviates from the corresponding result 
obtained in [271 ] . 

Authors of Ref. [3l[ calculated the susceptibility xv 
including only the zero-th component, xv.o, for vf = 0, 
performing an expansion for small q . They found the 
leading term oc q 2 /2 m. However, it has the opposite sign 
(see (48) in [3l[ ) to the second term, oc 1b , in (40) , that 
would yield reaction probability in case, if bare vertices 
were used. Note also that the key equations (35,38-45) in 
[3~L| differ from Larkin-Migdal equations (III) ( for T < 
A as supposed in [34| . and for vp = as assumed in 
[3~fl ] ) . As it follows from ([T4"]) and ([To]) our expression for 
(C(n, q, oc q 4 Vp n vanishes, if vp, n — > 0, although 

q 2 /m* terms were present in the original loop integrals. 



B. Emissivity on axial- vector current 

Now we focus on the process going on the axial-vector 
current. Then 



^5 



q 2 (vXA,i(n, q))n + (3w 2 -2q z ) (xA,o(n, q)). n 



u{<tX*A,i(n,q))n - uj ({qv)xA,o(n, q)) n 



(33) 



The last two crossing terms in the squared brackets can- 
not be eliminated. Keeping only terms oc vp we cast ([33]) 
as 



q 2 V F (L(n, q; +1)) H + (3 uj 2 - 2 q z )(L(n, q; -l)) n 



- u 2 (L(n,q;-l))fi - uj ((qv ) L(n, q; -l))n 

,2 



-a 2 pvpf [l + (l-|il)-|] % (^). 



(34) 



As in case with the vector current, simplifying we could 
put j A = 1 since -f A ~ 1 + 0(g% n v F ). 
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The contribution of the axial- vector current to the neu- 
trino emissivity is determined by 



(i + ^-^)4„Q (0n) H- 



(35) 



Second term in round brackets of f|35[) has been men- 
tioned already in [r| and then recovered in [HJ . Our 
coefficient (11/21) is twice larger than that presented 
in those works. We notice that the integral I s /2 = 
{u'v — v'u) 2 , where u and v are coefficients of the Bo- 
golyubov transformation, is in [HJ twice as large as that 
in [lj|. In agreement with the former evaluation, we ar- 
rive at the coefficient 11/21 rather than at 11/42, as pre- 
sented in [3, S3- The first term in (|35| (for m* = m) is 
the same, as in Ref. 22j, which calculated this relativistic 
correction for the fist time. The factor (m*/m) 2 does not 
arise in our calculations, since the mass renormalization 
is performed everywhere, including the vertices. Other- 
wise the Ward identity would not hold for the renormal- 
ized "bare" vertex rf, . The third term related to the 
time-space component product was not considered be- 
fore. 

Finally for the neutron PBF emissivity on the axial- 
vector current we obtain (for one neutrino flavor) 



iPBF 



* 2 2 
9a "F, n £ 



(On) 



(36) 



The resulting emissivity is the sum of contributions 
and (1361. 



11PBF _ pPBF 



nPBF 
t vu,A 



c iiPBF 
-vv,A ■ 



(37) 



The axial- vector term, being oc Dp, is now the dominating 
contribution. Thus, the ratio of the emissivity of the 
neutron PFB obtained here to the emissivity calculated 
in (l6Ll22l|. where main contribution was due to the vector 
current, is 



i?(nPFB) 



.nPBF 



(On) 



*2 2 
9A v F.n 



= F n V 



F,n ■ 



For n = no — 0.17 fm , m* = 0.8 to, we estimate 
F n ~ 0.9-1.2, v F . n ~ 0.36 and R ~ 0.12—0.15. For 
n = 2no, m* = 0.7m, R increases up to 0.24 — 0.32. This 
is in drastic contrast with estimations i?(nPFB) ~ 10~ 3 
in [13] (being actually valid only for the rate of par- 
tial vector current contributions, rather than for the full 
emissivities) and i?(nPFB) ~ 5 • 10~ 3 obtained in [311 ]. 



V. NEUTRINO EMISSION VIA PROTON PBF 



Now we turn to the proton PBF processes. If protons 
were the only particles in the neutron star medium, we 

could use the results obtained above for the neutron PBF 

» _ „{p) 
fly 



and just replace g v 

fpp and gnn * 9ppj &a 



C V, v F,n — ¥ Vp^p, f n , 



A. Emissivity on vector current 



Since the bare vertex yields now (fly'') 



0.002 



compared to (fl(T') 2 = 1 in the neutron case one could 
naively think that the emissivity of the proton PBF 
process on the vector current is suppressed by a factor 
~ 10" 3 (A p /A„) 13 / 2 e 2(A J'" A " )/T compared to the emis- 
sivity of the neutron PBF process. For imaginary purely 
proton matter we would find in the vector channel 



-w.y 



(Op) 2 J* 4 



4 

15^ 3 { T h 



(38) 
(39) 



In addition to that the emissivity is already suppressed 
in the vector channel by factor -^Vp , Leinson and 
Perez [27J found extra suppression. They included inter- 
action of protons via photons. This produces new con- 
tributions to the susceptibility x °f the following type 



Dotts assume infinite summation of the bubble chains 
with all four types of the vertices shown in The wavy 
line is the dressed photon Green function. Simplifying, 
Ref. [13] used the static Coulomb potential instead. To 
illucidate the origin of differences in our estimations with 
those in (27| we will use the same approximations. Effec- 
tively the summation leads to the replacement (2?| 



t~v.o 



Ty.o 
ecO?)' 



(40) 



where ec(q) — ^- + UJ pi/ uj2 is the dielectric constant, 
An e 2 n p /m* is the proton plasma frequency with e 2 



1/137. Setting m* = m, uj ~ 2A p , (w ~ 2 A + 0(T) for 
the PFB processes), A ~ 1.76T C and T c ~ 1 MeV for 
x. p = iip/n ~ 0.03 at n = n , cf. Fig. 2 of Q, we obtain 
ec (fl) — 1 + 0.3 ~ 1 in disagreement with estimation 
of [271 ] . where applying their result to the neutron star 
matter authors put n p = hq and u> ~ T c that resulted in 
the estimation ec(q) ~ 10 2 . Thus, a suppression factor < 
10~ 6 of the emissivity of the proton PBF process quoted 
in [13] is misleading. Note that correction of the vertex 
PO")) affects only the process on the vector current, since 
photon is the vector particle. 

For the neutron star matter, the replacement (|4T)|) is 
not sufficient, since protons are embedded into the elec- 
tron liquid of the same concentration and into a much 
more dense neutron liquid. Renormalization of the weak 
vector interaction of protons in this case can be taken into 
account, if we replace the bare coupling Ty as follows 



1 /j • 



(41) 
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Dots stand for other graphs not shown explicitly, like 
A(1232)n-loop term, etc. Simplifying, we ignore these 
rather small correction terms. The second graph has been 
incorporated in 0, [TtJ that results in the shift 

cv -*c v - /^pP _1 (no)5Ri„„7(/^„), 

where 7 _1 (/nJ = 1 - fnnP^ 1 { n oW L nn and L nn = 
{L(n,q;g = 0))^/a 2 = p(qv/(uj — qv))n is the Lindhard 
function. This correction (although oc Vp) leads to the 
strong enhancement of the tiny bare vertex. A numeri- 
cally larger correction comes from the electron-electron- 
hole polarization term (the third graph) . Such a possibil- 
ity has been discussed in 19] for the process of a possible 
massive photon decay and then it was taken into account 
in the proton PBF emissivity in [20']. Altogether these 
corrections can be incorporated to the resulting expres- 
sion for the emissivity of the proton PBF process with 
the help of the replacement, cf. [3], 



4 



T 



f% p p- 1 (n o nL n n 1 (f? m )+0.8C v 



(.42) 



Here C ve = 1 is the electron weak vector coupling. Thus 
we find 



pPBF _ (Op) T * 4 



(43) 



where the pre- factor T v ~ 1. Finally we obtain an esti- 
mate 



£ pPBF 

Rv\p/n] = - r 



4/3 ( £X P\ 

p uj 



A„ \ 13/2 



,2{A n -A p )/T 



(44) 



The ratio Rv[p/n] is sensitive to the values of the proton 
and neutron gaps as functions of density, A„. p (n), the 
proton fraction x p , and the temperature T. 



B. Emissivity on axial- vector current 

Now we consider axial-vector channel. Photon ex- 
change does not contribute in this channel. The main 
correction to the vertex comes from the iteration of the 
rm-loops: 



9pn 



(45) 



+ . . . 



Simplifying, we will suppress correlation factors, such as 

l 2 (9nn) - 1 + 0(<C<J- ThuS We ° btam 



c pPBF 
~vv,A 



(Op) ^ *2 2 
e vv rj 9 A v F,p ■ 



(46) 



Comparison with [221 ] can be done quite similar to that 
performed above for neutrons. 



We conclude 

pPBF pPBF pPBF ^ pPBF 

For the ratio R[p/n] we find 



(47) 



R\p/n] = 



pPBF e 



pPBF 
vv,A 



PBF — ,nPBF 
~-vv ^vu,A 



(A„-A p )/T 



(48) 



The ratio R[p/n] is sensitive to the choice of pairing gaps, 
temperature and the proton fraction x p and can be both 
^ 1 and > 1 . 



VI. CONCLUSIONS 

In this paper we re-calculated neutrino emissivity via 
neutron and proton pair breaking and formation pro- 
cesses. We used Larkin-Migdal-Leggett Fermi-liquid ap- 
proach to strongly interacting systems with the pairing. 
Compared to the Nambu-Gorkov formalism, the Larkin- 
Migdal-Leggett approach allows for different interactions 
in the particle-particle and the particle-hole channels, as 
it is the case for nuclear matter. 

To be specific we focused our discussion on the ISq 
pairing. We support statement of (2?| that medium ef- 
fects essentially modify vector current vertices. Only the 
careful account of these effects allows to fulfill the Ward 
identity and to protect conservation of the vector current. 
Compared to the emissivity calculated in [l6[ the partial 
contribution to the emissivity on the neutron vector cur- 
rent proved to be dramatically suppressed, roughly by 



a factor ~ 0.1 



x v 



F,n> 



is the neutron Fermi veloc- 



ity, cf. [2JJ. A similar suppression factor arises for the 
partial contribution to the emissivity on the proton vec- 
tor current, provided one replaces neutron Fermi velocity 
by the proton Fermi velocity. Electron-electron-hole and 
neutron-neutron-hole polarization effects play a crucial 
role in the latter estimation. Proton-proton-hole polar- 
ization effects are suppressed (these statements are at 
variance with the estimations in [271] ). 

The dominating contribution to the neutron and pro- 
ton pair breaking and formation emissivity comes from 
the weak axial- vector current. Finally, the neutron pair 
breaking and formation emissivity proves to be sup- 
pressed compared to that of [l(| by a factor ~ 0.12-0.15 
at nuclear saturation density and ~ 0.24-0.32 at twice 
nuclear saturation density. For the proton pair break- 
ing and formation the emissivity deviates only by a fac- 
tor close to unity from the expression used previo usly 
in [22j • These our findings differ from those in [27], l3lj |. 
where authors concluded that the neutron and proton 
pair breaking and formation emissivities are dramatically 
suppressed. Modifications of the neutron and proton pair 
breaking and formation reaction rates that we have found 
are probably not strong enough to essentially influence 
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on the values of surface temperatures of neutron stars 
computed previously. 

One may rise the question how much emissivity of 
other relevant neutrino processes might be changed, if 
medium effects in presence of nucleon pairing are cor- 
rectly included. Although we did not perform corre- 
sponding cumbersome calculations let formulate our con- 
jectures: 

In the reactions with charged currents, as the direct 
Urea and the modified Urea processes, the transferred 
neutrino energy is uj ~ /x e = pf, p ^> 2A. Therefore 
the anomalous Green functions are taken in the limit 
u> ^> 2A. In this limit ^-function tends to zero (as it fol- 
lows from the corresponding asymptotic in (fT5)) ). Effects 
of normal correlations and pion softening were evaluated 
in 0, IH EH E3 > resulting in significant enhancement of 
two-nucleon reaction rates. Specifics of the superfluid 
matter manifest themselves in reactions with charged 
currents mainly through the phase-space suppression fac- 
tors. 

The two-nucleon bremstrahlung processes going on the 
neutral currents are similar to the pair breaking and for- 
mation reactions. In a normal phase the emissivities are 
governed by the axial- vector current Thereby, we 

do not expect a strong suppression of these rates in a 
superfluid phase except the natural phase-space suppres- 
sion estimated in previous works. As in the case of the 
two-nucleon reactions going on charged currents, nucleon 
short-range correlations and pion softening significantly 
influence the reaction rates, cf. 0, [ll| . 

Our findings are relevant also for calculations 
of the quark-pair breaking and formation processes 
and other quark propagation processes in the color- 
superconducting medium, which use bare- loop results, 
e.g. see Refs. [H, Note that since the pairing 

gaps in the color superconductors can be rather large, 
2 A > /i e , both reaction rates on neutral and charged cur- 
rents (Urea) might be affected. Ref. [2(| considered neu- 
trino scattering off breaking pairs in color-flavor-locked 
medium within the Nambu-Gorkov formalism. However, 
they included only the zero-th component of the ver- 
tex, and their expressions for the current-current cor- 
relator in the non-relativistic limit do not coincide with 
the Larkin-Migdal-Leggett expressions that we have re- 
produced above. 

An interesting observation was made recently in 
Ref. [4f| . A natural explanation of the super-burst ig- 
nition would require a strong suppression of the neutron 
pair breaking and formation emissivity for low baryon 
densities. For n ~ 10 12 g/cm 3 we estimate a suppression 
factor as - 0.1 x (n/n ) 2/3 ~ 0.003. 

Another relevant issue is related to absorption and 
scattering processes of low-energy (uj ^ few MeV) neu- 
trinos and antineutrinos on nuclei. In absence of the 
electron Fermi sea the correlation effects may manifest 
themselves in reactions on both neutral and charged cur- 
rents. There are different sources for neutrinos of such 
energies, e.g., reactor neutrinos are good candidates. The 



Sun and supernova neutrinos also have pronounced low 
energy tails. The observation of supernova neutrinos 
might provide us with unique information on the core- 
collapse and on the compact star formation and cooling 
[46| . Geo-neutrinos and antineutrinos from the proge- 
nies of U, Th and 40 K decays inside the Earth bring to 
the surface information from the whole planet, about its 
content of radioactive elements [47l |. Finally, verifying 
the existence of the relic neutrino sea with temperature 
T„/T 7 = (4/11) 1 / 3 represents one of the main challenges 
of the modern cosmology [48[ . 

From a general point of view, our results strongly sup- 
port the conclusion of Refs. @, |3, EH 53 about the 
essential role played by different medium effects in the 
neutrino evolution of neutron stars, as it was demon- 
strated in the framework of the " nuclear- medium cooling 
scenario" [ill EH S3- Without a proper inclusion of 
medium effects it is difficult to reach sound conclusions. 
Further investigations in this direction are required. 
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APPENDIX A: CALCULATION WITH BARE 
VERTICES IN THE LOOP 

Let us now comment on results of previous calcula- 
tions of the neutron PBF emissivity. If, as in calcula- 
tions [ill H3 , one neglects contributions from the anoma- 
lous vertices r, i.e. M ^ and C — ► L in (fl"8]) , and puts 
7 a = 1 , one may recover results of calculations with bare 
vertices ("b.v."). Expression (|2"5)) for kv becomes 



„b.v. 



is 



q (L(n,q;+l))n 
+ v 2 ((lo 2 + (qn) 2 - q 2 ) L(n, q; -l)) n 
- to{(qv)(L(H,q;+l) + L(H,q;-l)))f^ .(Al) 

In the axial-vector channel the corresponding quantity is 
k\ v - cx Vp. Therefore we will keep also the v F corrections 
in (|A1[) . In calculations [rIUS] such terms were dropped. 
The second term in (|A1|) can be indeed neglected, as 
being cx v F , since (L(n, q; — 1))« oc v 2 ,. The third crossing 
term of the order v 2 ? disregarded in previous calculations 
can be dropped only, if conditions (f2"0"l) hold, that is not 
the case for the bare vector current-current correlator. 
We keep this term here. For \q \ up <C to (i.e. for vp <C 1 
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since \q \ Z; uj) we get 

+ (A2) 

The first line in (|A2|) comes from the expansion of the first 
term in (|Alj) . The first term in the second line follows 
from the crossing term in (|A1[) . Using (|A2|) we calculate 




Dropping the Up corrections and using expression Ij27p for 

Qy^ one may reproduce the vector current contribution 
to the emissivity obtained in Refs. [16|, |22j (and [2L 1 1 71] . 
provided one sets there jv = 0). 



The expression for the emissivity of neutron PBF on 
the axial- vector current calculated with the bare vertices 
in the loop does not deviate from that in (|35p derived 
above with the full vertices. Note, however, that second 
and third terms in (|35[) differ from those used in previous 
calculations . 
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